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Abstract 



LOPES is a digital antenna array at the Karlsruhe Institute of Tech- 
nology, Germany, for cosmic-ray air-shower measurements. Triggered 
by the co-located KASCADE-Grande air-shower array, LOPES detects 
the radio emission of air showers via digital radio interferometry. We 
summarize the status of LOPES and recent results. In particular, we 
present an update on the reconstruction of the primary-particle prop- 
erties based on almost 500 events above 100 PeV. With LOPES, the 
arrival direction can be reconstructed with a precision of at least 0.65°, 
and the energy with a precision of at least 20 %, which, however, does 
not include systematic uncertainties on the absolute energy scale. For 
many particle and astrophysics questions the reconstruction of the at- 
mospheric depth of the shower maximum, X max , is important, since it 
yields information on the type of the primary particle and its interac- 
tion with the atmosphere. Recently, we found experimental evidence 
that the slope of the radio lateral distribution is indeed sensitive to 
the longitudinal development of the air shower, but unfortunately, the 
A" m ax precision at LOPES is limited by the high level of anthropogenic 
radio background. Nevertheless, the developed methods can be trans- 
ferred to next generation experiments with lower background, which 
should provide an X mlLX precision competitive to other detection tech- 
nologies. 

1 Introduction 

Radio emission from air showers is studied both experimentally and theo- 
retically already for about 50 years, and in the last years significant progress 
has been achieved. Analyses with digital antenna arrays like LOPES could 
show that at least the energy and the arrival direction can be reconstructed 
precisely enough to make radio arrays a useful tool for cosmic-ray physics 
above 10 eV. These results are generally confirmed by other experiments, 
e.g., CODALEMA [H [2] and AERA [3J gj. It still remains to be demon- 
strated that radio measurements can contribute significantly to disentangle 
the composition of ultra-high energy cosmic rays. Recently, we could show 
that LOPES measurements are sensitive to the longitudinal shower devel- 
opment, and thus - at least on a statistical level - to the composition of 
the primary cosmic rays. However, the currently achieved precision on the 
atmospheric depth of the shower maximum, X majx , is insufficient due to the 
high level of anthropogenic radio background at the LOPES site. 

Since the LOPES experiment with its different setups has been described 
in several previous publications (e.g., reference [5j [H [7]), we will give only 
a short summary of our setup and the analysis procedures, here. LOPES 
is a digital antenna array co-located with and triggered by the KASCADE- 
Grande experiment [8] at the Karlsruhe Institute of Technology in Germany. 
It was built in 2003 as LOFAR [9] prototype station and made the proof-of- 
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principle that air showers can be measured with digital radio interferometry 
0. 

Because of its success, LOPES was enhanced several times. In 2005, it 
was extended from 10 to 30 inverted v-shape dipole antennas, which were 
all aligned in the east-west direction, because we assumed that the radio 
signal is mainly east-west polarized. End of 2006, half of the antennas 
were rotated to north-south alignment, and thus we could test the previous 
assumption - it is true at least for most shower geometries [10]. In 2009, we 
started LOPES-3D [7], which consists of 10 tripole antennas (= three crossed 
dipoles in the same location), to study the polarization in more detail. 

Generally the analyses and developments performed at LOPES aim at 
three main goals: 

Reconstruction of air shower parameters: 

Already in the 1960s historic experiments based on analog techniques demon- 
strated that in principle the arrival direction and the primary energy can 
be reconstructed from radio measurements [11 j . but they could not demon- 
strate a precision competitive to other air-shower detection techniques. This 
now has changed, and we will give upper limits for the precision in the next 
section. In addition, we have demonstrated that a reconstruction of X max 
is possible |12j . and current analyses aim at improving the reconstruction 
methods. 

Understanding of the radio signal: 

To maximize the benefit of radio measurements for air-shower and cosmic- 
ray physics, a sufficient understanding of the radio emission process is re- 
quired. LOPES [5] and CODALEMA [13] confirmed that the dominant 
emission mechanism is the geomagnetic deflection of air-shower particles. 
Currently, we try to further improve the understanding of the radio emis- 
sion in two ways: First, we analyze measured properties of the radio signal, 
like the polarization [T5], the lateral distribution [15] . and the wavefront 
shape [16J . Second, we compare measurements to Monte Carlo simulations 
of the radio emission. This way we found strong indications that on top of 
the geomagnetic deflection and the sub-dominant Askaryan effect |17j (= 
net charge variation), also the refractive index of the air has an important 
influence on the radio signal, since it changes the coherence conditions for 
the emission [18j Qjj] . 

Prototyping for large scale radio arrays: 

Aside from the physics goals, LOPES aims at technical developments for 
future, large scale radio arrays. Different antenna types [20], self-trigger 
algorithms, analog and digital radio electronics have been developed and 
tested within the LOPES STAR extension [21] . The results have influenced 
current state-of-the-art experiments like AERA [3] and Tunka-Rex [22], and 
now these efforts are stopped at LOPES and continued in these experiments. 
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Figure 1: Shower cores of events measured with the east-west aligned LOPES 
antennas within the fiducial areas of the KASCADE particle-detector array 
respectively its high energy extension 'Grande'. 



In this paper we focus on the first goal, i.e. the reconstruction of air- 
shower parameters. The progress in other topics is covered by LOPES pro- 
ceedings in this issue [TU [TU], or was already published in journals (see 
www . lopes-pro j ect . org for an overview on LOPES publications). 



2 Results 

Since the first studies on the LOPES precision for the arrival direction [23J 
and the energy [23] we made significant progress by: a) improving our am- 
plitude and time calibration [23L2S], b) introducing systematic, quantifiable 
quality cuts for the event selection [27], c) improving many details in our 
analysis software, e.g., up-sampling, and a better fit algorithm for the direc- 
tion reconstruction. Due to the considerably higher statistics we used the 
signal of the east-west aligned antennas, only, which is analyzed under the 
approximation that the radio signal is purely east-west polarized. We are 
planning to study whether we can improve the reconstruction precision by 
including measurements of the differently aligned antennas. Nevertheless, 
already with the present analysis we can give convincing upper limits for the 
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direction and energy precision, and are able to reconstruct the atmospheric 
depth of the shower maximum, X max . 

The reconstruction of shower parameters is based either directly on cross- 
correlation (CC) beamforming, or on the amplitude and arrival time in each 
individual antenna - after identifying the pulse with the CC-beam. During 
CC-beamforming the traces recorded by different antennas are shifted in 
time and combined to one single quantity by calculating a cross-correlation. 
To calculate the expected arrival times in each antenna for a given direc- 
tion we assume either a spherical or a conical wavefront. The CC-beam is 
maximum when it is formed into the direction of the radio signal, i.e. when 
the time shift corresponds to the arrival time difference of the radio signal 
in the individual antennas (for more details see references [28|, IB"]). 

In most cases the measurement precision is limited by noise [29J, and 
only at high signal-to-noise ratios by the amplitude and time calibration of 
LOPES. For high signal-to-noise ratios, the uncertainty of individual am- 
plitudes is 5 % when comparing different antennas, and another 5 % when 
comparing different events with each other, due to the effect of changing 
environmental conditions. In addition there is a scale uncertainty on the 
absolute LOPES amplitude of about 35%, which is relevant when compar- 
ing LOPES measurements to simulations or other experiments. The rela- 
tive timing accuracy is in the order of 1 ns and in almost all cases negligible 
against the noise induced uncertainty on the pulse time of up-to 18ns [29] . 

For the following analyses we selected events whose energy was recon- 
structed above 100 PeV by KASCADE, respectively KASCADE-Grande. 
Furthermore, we require a clear radio signal in the CC-beam, and several 
standard quality cuts [27]: e.g., the core must be inside of the fiducial areas 
of KASCADE (281 events), respectively KASCADE-Grande (204 events). 
Figure Q] displays a map of the cores of the selected LOPES events. 

2.1 Direction 

For the direction reconstruction we used the KASCADE(-Grande) direction 
as input, and then searched within a 2.5° neighborhood for the maximum of 
the CC-beam. We made several cross-checks that the 2.5° neighborhood is 
sufficiently large to find the global maximum of the CC-beam: First, 2.5° is 
large against the KASCADE-Grande angular resolution (<< 1°). Second, 
we tried larger neighborhoods and, third, randomly varied the input direc- 
tion by a value about three times larger than the LOPES angular resolution, 
but the results did not change significantly. 

The total angular difference between the reconstructed LOPES direc- 
tion and the KASCADE(-Grande) input direction is given in figure [2j We 
take the mean difference as upper limit for the LOPES direction precision 
for cosmic-ray air showers. It is 0.72 ± 0.02° for the spherical wavefront 
assumption and 0.65 ± 0.02° for the conical wavefront assumption. In addi- 
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Figure 2: Deviation of the KASCADE, respectively the KASCADE-Grande 
direction and the LOPES direction reconstructed with a spherical and a 
conical wavefront assumption. The mean and its statistical uncertainty has 
been determined by fitting a Gaussian. 



tion to the analysis we presented in reference [16] , this is a further indication 
that the 'true' radio wavefront is better approximated by a cone than by a 
sphere. 

2.2 Energy 

The energy of the primary particle can be reconstructed from the amplitude 
of the radio signal in several ways. By comparing the LOPES energy recon- 
struction with the KASCADE- (Grande) result we found that at least the 
statistical precision is similar or better than the KASCADE-(Grande) preci- 
sion. However, we cannot determine the LOPES precision accurately since 
the KASCADE- (Grande) precision is not known accurately for individual 
events, and because we have a selection bias. Therefore, we determine only 
upper limits on the LOPES energy precision. 

For the LOPES energy reconstruction we use different methods which 
all are essentially modifications of the formula presented by Allan 

I ^1 

-^Allan OC - --— — - (1) 

sin a ■ cos a ■ exp (— a^is/do) 

with E the energy, e the radio amplitude, a the geomagnetic angle (= angle 
between the shower axis and the Earth's magnetic field), 9 the zenith an- 
gle, riaxis the distance to the shower axis, and do a slope parameter for an 
exponential lateral distribution. 

One method relies on the amplitude measurements in the individual 
antennas, and determines the amplitude at a typical distance by fitting the 
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Figure 3: Comparision of the LOPES reconstructed energy with the 
KASCADE, respectively the KASCADE-Grande energy. The dashed lines 
indicate an energy shift by ±20 % which roughly is the KASCADE(-Grande) 
energy precision. 



lateral distribution [30] (see reference [31] for new LOPES results based on 
this method). Another method already used for earlier LOPES results |24J 
is based on the amplitude of the CC-beam, i.e. one single quantity per event 
combining the measurements of all antennas. We present and update here 
(figure [3]), using the following formula for energy reconstruction: 

„ . e C c,Ew/31MHz 

£lopes,cc = const ■ — (2) 

\v x B\ew ■ exp (-d axis /180m) 

with ecc the amplitude of the cross-correlation beam which we divide by 
the effective bandwidth of LOPES (43 - 74 MHz), and \v x £?| E \y<i the 
east-west component of the geomagnetic Lorentz force unity vector, which 
is different from sin a in most cases. Unlike Allan, we did not divide by 
cos 9, since this would increase the spread in figure [3l and thus would make 
the energy reconstruction worse. The proportionality constant has been 
determined by a cross-calibration with KASCADE(-Grande) such that the 
mean deviation between the LOPES and KASCADE(-Grande) energy is 
0. The results are constxASACDE = 17PeV-m-MHz/uV, and constcrande = 
13 PeV-m-MHz/uV. There are two possible explanations for the difference 
between both constants: First, in contrast to the Grande energy reconstruc- 
tion, the KASCADE reconstruction has not been designed for this energy 
range. This might also explain the non-linearity at high energies in figure 
[3] (left panel). Second, the Grande and KASCADE events are measured at 
different axis distances (cf. figure [1]), and the correction by an exponential 
lateral distribution with a fixed scale factor might be oversimplified to cover 
the full distance range. 
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Figure 4: Relative deviation of the LOPES reconstructed energy and the 
KASCADE, respectively the KASCADE-Grande energy. The width and its 
statistical uncertainty has been determined by fitting a Gaussian. 



Nevertheless, the used energy reconstruction formula is sufficiently good 
to state an upper limit for the LOPES energy precision of 20% (figure 
U]), at least for the KASCADE events, i.e., events whose core typically is 
inside the LOPES array. Of course the total uncertainty is worse, since 
systematic uncertainties on the energy scale have to be included (see Ref. [32] 
for a discussion on the systematic uncertainties of the KASCADE-Grande 
energy reconstruction we used for cross-calibration). In addition, the reason 
for single outliers like the two events in figure [3] (right panel) has to be 
understood: Their geomagnetic angle is not exceptionally small, but they 
are measured at a relatively large distance from the LOPES array of several 
100 m, and the used energy reconstruction might not be reliable in this case. 

2.3 Shower Maximum 

Recently, we provided an experimental proof that radio measurements are 
indeed sensitive to the longitudinal shower development, as has been the- 
oretically assumed since long. We observe a correlation between the mean 
pseudorapidity of high energy muons measured by the KASCADE muon 
tracking detector [33J and the slope of radio lateral distributions measured 
by LOPES [12]: young showers typically have a large mean muon pseudo- 
rapidity and a steep lateral slope, and old showers vice-versa. In reference 
[3T] we describe a method how to reconstruct the atmospheric depth of the 
shower maximum, Ajnax? via the lateral slope. Using REAS3 simulations 
|34j for calibration, we achieve an upper limit for the LOPES X max pre- 
cision of 90g/cm 2 , and the simulations indicate that in principle, i.e. in a 
situation with negligible noise, a precision of better than 30g/cm 2 should be 
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Table 1: LOPES reconstruction precision achieved with present analysis 
techniques (not including additional systematic scale uncertainties). 



Reconstruction quantity / method 


LOPES precision 


Direction: 




Spherical CC-beamforming 


< 0.72° 


Conical CC-beamforming 


< 0.65° 


Energy: 




Amplitude of CC-beam 


< 20% 


Amplitude at typical axis distance 


< 20% 


Amax- 

Slope of lateral distribution 


90g/cm 2 


Cone angle of radio wavefront 


ps 200g/cm 2 



possible. Thus, we conclude that the main reason for the large measurement 
uncertainty is the high ambient noise level at the LOPES site. 

We have tested an independent, second method for the reconstruction 
of X max . REAS3 simulations show that the cone angle of the wavefront 
is proportional to X max after a correction for the shower inclination. Ac- 
cording to the simulations, also for this method a precision of better than 
30g/cm 2 should be achievable, and also in this case the real precision of the 
LOPES measurements is much worse |16j . Nevertheless, a combination of 
both methods might improve the measurement precision such that we are 
able to distinguish at least heavy from light nuclei, even in the noise envi- 
ronment of LOPES. In principal, there is a third method we have not tried 
at LOPES, since also the frequency spectrum of the radio signal ought to be 
sensitive to the longitudinal shower development [35J. However, due to the 
large noise-related uncertainties of LOPES we could only demonstrate on a 
statistical basis that the radio amplitude decreases towards high frequencies 
[3SIEZ], but we are not able to determine the spectral slope for individual 
antennas on a per-event basis. 

3 Conclusion 

LOPES has been contributing to the development of the radio measurement 
technique for air showers for almost one decade, and still makes significant 
progress. We have a total data set of more than 500 high quality radio 
events which are used in many ways. Comparisons between new simulations 
and LOPES data reflect the improved understanding of the radio emission 
[19] . and comparisons between the LOPES and KASCADE-Grande recon- 
struction of shower parameters reflect the potential of the radio method in 
general. For the shower direction and the energy of the primary particle, 
LOPES demonstrates that the radio technique can compete with the preci- 
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sion of established techniques. For X max , we developed promising methods, 
but the actual precision is limited by the high level of human-made back- 
ground at the LOPES site. Nevertheless, experiments in regions with lower 
background like AERA [3J, LOFAR [9] or Tunka-Rex [22J should be able to 
achieve a precision competitive to air-fluorescence and air-Cherenkov mea- 
surements. These experiments can also test theoretical predictions that with 
radio arrays an energy precision below 10% should be achievable [50"1 13T]. 

Despite all this enthusiasm, our LOPES results indicate also that the 
stand-alone operation of radio arrays is difficult - at least in radio-loud en- 
vironments. Self-triggering on the radio signal turned out more complicated 
than previously thought, since air-shower pulses are hard to distinguish from 
anthropogenic interferences. In addition, the absolute scale of the radio am- 
plitude is still only understood up-to a factor of 2, which would affect both 
the energy scale uncertainty of stand-alone radio arrays as well as simula- 
tions of their efficiency. Nevertheless, digital radio antennas seem to be a 
promising extension for particle detector arrays, since radio arrays have a 
duty-cycle of almost 100%, and most likely can increase the precision of 
particle detector arrays for the energy and composition of ultra-high energy 
cosmic rays. 
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